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§ ! Abstract 

<N 

We present a theory for the generation of mesoscale (kr g <C 1, where r g is the cosmic ray 
gyroradius) magnetic fields during diffusive shock acceleration. The decay or modulational insta- 
bility of resonantly excited Alfven waves scattering off ambient density perturbations in the shock 
environment naturally generates larger scale fields. For a broad spectrum of perturbations, the 
physical mechanism of energy transfer is random refraction, represented by diffusion of Alfven 
wave packet in k— space. The scattering field can be produced directly by the decay instability or 
by the Drury instability, a hydrodynamic instability driven by the cosmic ray pressure gradient. 
This process is of interest to acceleration since it generates waves of longer wavelength, and so 
enables the confinement and acceleration of higher energy particles. This process also limits the 
^ ' intensity of resonantly generated turbulent magnetic field on r g scales. 
O " 

^: 

q . 1. Introduction and Overview 

^ ■ 1.1. Status of the CR Acceleration Problem 



There is an increasingly popular view that the CR spectrum above the "knee" (about 10 15 eV) 
is produced by the same accelerator (or accelerators) as the part below the knee. The latter portion 
of the spectrum is better understood. As is usually argued, spectrum softening at the knee (from 
E~ 2J to E~ 3A ) forces adherents to the premise of extragalactic origin for CRs above the knee to 
also explain why the galactic part of the spectrum terminates at precisely the point where the extra- 
galactic part appears. Therefore it appears desirable to explain the CR spectrum straight up to its 
next feature, the "ankle" (at about 10 18 eV) by a single mechanism operating in one accelerator or 
in a group of similar accelerators. SNR shocks are now considered to be the most promising can- 
didate for that purpose. Recent discussion of available suggestions that seek to accelerate cosmic 
rays (CRs) in supernova remnants (SNR) to energies beyond the "knee" can be found in Parizot 
et al. (2004); Hillas (2005). 

Of course, SNRs are almost certainly responsible for the CR (at least electron) acceleration 
below the knee, which is documented in several ways (see Aharonian, 2004 for a comprehensive 
review of detection techniques and physical processes). In particular, there is an evidence (Koyama 
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et al. 1995; Tanimori et al. 1998; Allen, Petre, & Gotthelf 2001; in the form of both synchrotron 
and inverse Compton radiation) that electrons of energies up to 1007W are accelerated in the 
supernova shock waves. Recently, accelerated electrons with record energies up to at least 450TeV 
have been discovered in a young shell-type supernova remnant G12. 8-0.0 (Brogan et al. 2005). 
Interestingly enough, this source was previously detected in the y-ray band by the HESS Cerenkov 
telescope (Aharonian et al. 2005). 

By association, the electron acceleration mechanism should be also responsible for acceler- 
ating the main CR component, namely the protons. Indeed, at ultrarelativistic energies particle 
dynamics is determined by particle momentum (i.e., rigidity) regardless of mass. However, the 
acceleration of protons in supernova remnant (SNR) shocks has not been confirmed conclusively. 
The only SNR where a signature of accelerated protons was claimed to be observed is the RXJ 
1713.7-3946 (Enomoto et al. 2002). This source was first detected in the y- ray band by the 
CANGAROO team (Muraishi et al. 2000), and later confirmed by HESS (Aharonian et al. 2004). 
However, the claim made by the CANGAROO team (Enomoto et al. 2002) became a subject of 
significant controversy (Reimer & Pohl 2002; Butt et al. 2002), primarily on the grounds that ex- 
planation of the data would require a break in the energy spectrum of accelerated protons, while 
the standard diffusive shock acceleration (DSA) mechanism does not predict the necessary break. 
A possible resolution of this controversy, based on nonlinear effects in the diffusive acceleration 
mechanism, along with the generation of a break due to modification of particle confinement in 
molecular clouds surrounding the remnants (via neutral particle effect on Alfven waves) was re- 
cently suggested by Malkov, Diamond, & Sagdeev (2005). The latest high resolution observations 
with HESS (Aharonian et al. 2006) are also consistent with the presence of a break in the spectrum 
in the TeV energy range, but only future, more accurate data will be able to distinguish between 
the different possible functional forms of the high energy spectrum predicted by different mod- 
els. The simplest (and by far the most popular) approach, based on a power-law spectrum with 
an exponential cut-off, produces almost as good a fit to the current data as the spectrum with a 
break (Aharonian et al. 2006) does. However, the spectral break is a distinguishing feature, that 
reveals valueable information about the physical processes responsible for its formation and about 
acceleration in general. These processes include, but are not limited to, quasi-abrupt changes in 

1. the dynamics of waves that confine particles (Malkov, Diamond, & Jones, 2002; Diamond 
& Malkov, 2004; Malkov et al., 2005) 

2. the overall acceleration regime (Drury, van der Swaluw, & Carroll 2003), and 

3. the particle confinement regime (Malkov & Diamond 2006). 

Further broad band studies of RXJ 1713.7-3946 are also needed to fully resolve the dilemma 
of hadronic versus leptonic origin of the TeV emission. Recent TeV observations with HESS 
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(Aharonian et al. 2004) allowed authors to suggest that both leptonic and hadronic components 
are being accelerated in this object. Furthermore, the latest observations with the same instrument 
(Aharonian et al. 2006) seem to indicate a decline in the proton spectrum at energies ~ lOOTeV if 
the observed spectrum can be interpreted as a by-product of the interaction of accelerated protons 
with the ambient gas. From a theoretical standpoint, a very interesting aspect of this spectrum is 
that it is significantly softer than one would expect from a strong shock in a nonlinear (spectral 
index~ 1.5) or even in a linear regime (spectral index= 2.0). Should future measurements confirm 
that it must be interpreted as a piece of a broken power-law spectrum (Malkov et al. 2005), rather 
than its cut-off (Berezhko & Voelk 2006), a new interpretation requiring new DSA physics, along 
the lines of (1-3) above will be necessary. Note that current numerical models (e.g., Berezhko & 
Voelk, 2006) do not include the wave-particle interaction effects necessarily involved in (1) and 
(3). Unfortunately, even the high quality TeV data of (Aharonian et al. 2006) degrades statistically 
at photon energies above 10 TeV, where the spectrum begins to decline. Nevertheless, the spectrum 
appears more like a broken power-law at 10 TeV rather than a power-law with an exponential cut- 
off. The latter is more consistent with the time-dependent acceleration or energy loss scenario (see, 
e.g., Fig. 12 in Aharonian et al., 2006, and Fig. 8 in Berezhko & Voelk, 2006) implemented in most 
numerical models. 

Interestingly, even observations of the background CRs pose similar problems to the standard 
acceleration theory. The most notorious one is perhaps the overall CR spectrum itself which, 
as discussed earlier, is too steep in its high energy part (i.e., above the "knee") to be explained 
straightforwardly by standard acceleration theory. Leaving aside energy dependent propagation of 
CRs as well as their reacceleration (which may also significantly influence their spectrum, as noted 
by Ptuskin et al., 1997), one can distinguish two problems. The first problem is the physical origin 
of the break ("knee") and the second problem is the very stringent requirements on the acceleration 
parameters, such as the turbulent component of the magnetic field, in order to reach to the very high 
"ankle" energy of ~ 10 18 eV. 



1.2. Approaches to Enhanced DSA 

The important quantities that regulate diffusive shock acceleration (DSA) are the strength 
and spectral distribution of the turbulent magnetic field, 85. The magnetic turbulence confines 
accelerated particles to the shock front by pitch-angle scattering and is believed to be produced 
by the particles themselves, via CR-Alfven wave resonance when CRs stream ahead of the shock. 
Within standard quasi-linear theory (which is strictly valid only for 85 <C Bo), pitch angle scattering 
of relativistic particles by the waves proceeds at the rate 
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where Q. and p are the (nonrelativistic) gyrofrequency and momentum. Resonance requires kr g = 
const . Particle self-confinement along the field is diffusive, and the diffusivity is K ~ c 2 /v. The ac- 
celeration time scale can be estimated as x acc ~ k/w 2 , where u s is the shock speed. The fluctuating 
part of the field is usually assumed to saturate at the level of the ambient field, 85 ~ Bq, which thus 
produces pitch-angle scattering at the rate of Q. (the gyrofrequency), which limits the particle mean 
free path (m.f.p.) along the field to a distance of the order of gyroradius. This constitutes the so 
called "Bohm diffusion limit". Under these circumstances, the mean field Bq sets the acceleration 
rate and the maximum particle energy. The latter can be expressed through the work done by the 
induced electric field (u s /c)Bq on the particles while they are carried with the shock at the speed 
u s over the length scale of the shock radius, R s . Thus, the maximum energy is E max ~ (e/c)u s BoR s . 
Thus, one arrives at the maximum CR energy accelerated in a typical SNR shock of about 10 15 eV, 
which is close to the "knee" but is three orders of magnitude below the "ankle". Before discussing 
any approaches to enhanced, beyond-the-"knee", acceleration, it is important to note that, due to 
the resonance condition kp = const, confinement of higher energy particles requires that longer 
waves must be excited. Put another way, inverse scattering or transfer of Alfven wave energy 
excited at kr g ~ 1 to longer scales is clearly beneficial to confinement and acceleration. 

In order to reach the energy of the "ankle" several suggestions have been made. One approach 
is to invoke the generation of a fluctuating field component 85 significantly in excess of the un- 
perturbed field Bq (Lucek & Bell 2000). Physically, such generation is deemed possible since the 
free energy source is the pressure gradient of accelerated particles, which may reach a significant 
fraction of the shock ram energy. Specifically, a free energy limit on the wave energy density 
(85/5o) may be related to the partial pressure P c of CRs that resonantly drive the waves by the 
relation (MacKenzie & Voelk 1982) 



Here Ma = u s /Va 3> 1 is the Alfven Mach number and pu 2 s is the shock ram pressure. Of course, 
when 85/5o exceeds unity, particle dynamics, and thus the particle confinement and acceleration 
rates depart radically from the quasi-linear picture that underpins the usual DSA theory and mod- 
eling. The simple case of a monochromatic wave with arbitrary 85/5o, in which particle dynamics 
are exactly integrable (Lutomirski & Sudan 1966), provides an important clue to the general case 
{e.g., Malkov, 1998). A critical parameter is kr* where k is the wave number and r* is the particle 
Larmor radius calculated with the perturbed 85 ^> Bq field, rather than Bq, Fig.l. Particles with 
kr* < 1 perceive a strong local field which is perpendicular to Bo and therefore their confinement 
in the Bq direction is good. Particles with kr* > 1 perceive only the averaged, rather than the local 
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field, which is weak. The 85 component (even if it is large) exerts only a (rapidly) oscillating force 
on these particles, which thus can escape along Bq. Put more mathematically, for kr* S> 1, the 
particle response is the Boltzmann response. For the resonantly driven waves kr g ~ 1, and so if the 
waves grow nonlinearly until r* <C r g , the particles that initially destabilize the waves are trapped 
by the wave, thus saturating the instability in the wave band corresponding to their energy (Lucek 
& Bell 2000). This is analogous to the saturation of the beam-plasma instability for the so called 
hydrodynamic regime, for which the wave stops growing when its amplitude is sufficient to trap 
the beam particles. The numerical studies by Lucek & Bell (2000) showed that at least for the case 
of an MHD background plasma and rather narrow wave (and particle energy) band, the amplitude 
of the principal mode can reach a few times that of the background field. Moreover, Bell & Lucek 
(2001) argue that in the case of efficient acceleration, field amplification may be even stronger, 
reaching a mG level from the background of a few /uG ISM field, thus providing acceleration of 
protons up to 10 17 eV in SNRs. In view of the above discussion, however, we would like to em- 
phasize that the simulations of Bell and Lucek are intrinsically limited precisely because they are 
narrowband. Thus, wave-particle interaction is quite restricted, and the simulations do not allow 
any interaction of Alfven waves with other types of fluctuations which, in reality, are likely to be 
present in a shock environment. Moreover, Bell and Lucek did not address the question of how to 
compute the particle scattering rate for 85 3> Bq. Certainly, it is not correct to simply plug such a 
large amplitude into the usual quasi-linear diffusion coefficient, equation (1). 

Recently, Ptuskin & Zirakashvili (2003) approached this problem from a different perspec- 
tive. They considered a Kolmogorov-type turbulent cascade to small scales, assuming the waves 
are generated by efficiently accelerated particles on the long-wave part of the spectrum. The ques- 
tion of why such a Kolmogorov-like cascade model is relevant to a shock environment was not 
addressed. They obtained a maximum particle energy similar to the result of Bell & Lucek (2001). 
Yet another consequence of high magnetic field fluctuation levels, which is based on the change 
of the expansion regime of the SNR shocks, is discussed by Drury et al. (2003). An important 
question that still remains, is how realistic the saturation level given by equation (2) isl Earlier 
studies by Voelk, Drury, & McKenzie (1984) and Achterberg & Blandford (1986) suggest that 
due to particle trapping the instability saturates at levels 85 ~ Bq. For Alfven turbulence, consid- 
eration of random parallel scattering (i.e., turbulent mirroring) (Achterberg 1981; Achterberg & 
Blandford 1986) implies a similar saturation level. Thus, the questions of the saturation level and 
the mechanism of confinement of high energy particles remain unanswered. 

Motivated by the above issues and by the problems in describing the observed high-energy 
spectra as outlined in the Introduction, Malkov & Diamond (2006) suggested a faster-than-Bohm 
rate acceleration which is also intimately related to the knee and other spectral break phenomenon. 
The mechanism does not require 6B/Bq^> 1 magnetic field fluctuations, since particles gain energy 
by bouncing between the scatterers convected with the gradually converging upstream flow in a 



-6- 



nonlinearly modified shock precursor. This is different than the bouncing between the upstream 
and downstream scattering centers, which is usually assumed in DSA models. Note that in the 
latter case the acceleration rate decreases with energy, on account of increasing particle mean 
free path (particle diffusion length at the shock). This generic deficiency of the DSA is not relevant 
to the recently proposed mechanism, since the increase in m.f.p. is exactly compensated by the 
increased flow-induced compression of the scattering centers between which the particle bounces 
in the precursor flow. The maximum momentum is estimated to be 

c L 

Pmax~ -j-P* (3) 

where L/L p is the ratio of the distance between the scattering centers (weak shocks) to the precur- 
sor length (roughly estimated to be < 10) and is the maximum momentum achieved during the 
standard phase of the DSA, which becomes the break point or "knee" of the final spectrum. The 
spectral index between and p max is steeper than the "standard" p~ 4 and its slope depends on 
details of particle interaction with scatterers. The acceleration time is 



W- (Pmax) ~ TiVL (p* ) In (4) 

P* 

where x NL (p) ~ 4k^ (p) / u 2 sh (with the Bohm diffusion K s ) is the nonlinear acceleration time 
(Malkov & Drury 2001) which only slightly differs from the upstream contribution to the stan- 
dard linear acceleration time. 



1.3. Inverse Cascade and Enhanced DSA 

The theory of enhanced acceleration via bouncing between scatterers undergoing compression 
in the precursor flow described above does not specifically address magnetic field effects. Indeed, 
one is naturally motivated to ask whether it is, in fact, possible to achieve 85 3> Bq in a turbulent 
environment where many types of different wave interactions are possible. If such high levels 
are not achieved, one then must confront the issue of how to confine high energy particles to the 
shock. To this end, Diamond & Malkov (2004) previously suggested an acceleration scenario 
in which the magnetic field may interact strongly with the shock as a result of the acceleration 
itself, which, in turn, may in fact be strongly enhanced. The mechanism of such enhancement is 
based on the transfer of magnetic energy to longer scales via wave-wave interaction, which we call 
inverse cascade for short, even though specific mechanisms of such transfer may differ from what 
is usually understood as a local, self-similar cascade in MHD turbulence. This transfer is limited 
only by an outer scale L out - likely the shock precursor size K(p max )/u s ~ r g {p max )c/u s > r g (p max ). 
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The conceptual picture of this process is simple. There are two populations of fluctuations naturally 
native to a shock environment. These are 

1. Alfven waves, resonantly generated at small scales, i.e., kr g ~ 1 

2. acoustic waves and density fluctuations, at kL 0Ut > 1 but not 3> 1. Acoustic modes can be 
generated by many processes. The density perturbation field naturally refracts the Alfven 
wave field, i.e., 



For a random array of scatterers, k evolves diffusively, so diffusion of the Alfven wave spectrum 
results, with 



This diffusion naturally spreads the Alfven wave population in k, and so both prevents the devel- 
opment of a narrow spatial band or beam with SB/Bq ^> 1. This scattering also generates longer 
waves, which can in turn confine higher energy particles. Here, the Alfven wave population density 
flux is simply 



where N is the Alfven wave population density (i.e., wave action density N = e/cfl). Since dN/dk > 
for kr g < 1 (i.e., Alfven waves are excited at kr g < 1), the flux is toward larger scales and 
lower k's. It is important to note that this process is non-local in wavenumber, and thus not a true 
cascade in the traditional sense, and also is not a 'dynamo', but rather a re-distribution of magnetic 
energy among different scales. This approach is in distinct contrast to the models of Bell & Lucek 
(2001) and Ptuskin & Zirakashvili (2003) discussed above, which work with magnetic fields on 
scale lengths of order of Larmor radius r g (p max ) of the highest energy particles, and even smaller. 
The advantage of an 'inverse cascade' for acceleration is that longer waves confine higher energy 
particles and that the turbulent field at the outer scale $B(L 0Ut ) = B rms (generated by the flux in 
wave number space) which necessarily must have long autocorrelation time can likely be regarded 
as an "ambient field", so far as accelerated particles of all energies are concerned. If B rms 3> Bq, 
then the acceleration can be enhanced by a factor B rms /Bo. Note that the resonance field §B(r g ) 
remains smaller than B rms , so that standard arguments about Bohm diffusion apply. Since the 
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fluctuations around the new 'background' field B rms remain relatively weak, they are not likely to 
dissipate so rapidly via nonlinear processes, such as induced scattering on thermal protons, as it 
is to be expected in the case 85 > Bq when the transfer of energy from resonant fluctuations 85 
to the large scale B rms is not taken into account. In particular, such nonlinear processes were not 
included in the enhanced acceleration models of Lucek & Bell (2000); Bell & Lucek (2001). These 
can reduce the resonant fluctuation field significantly. 

As should be clear from the discussion above, an adequate description of the acceleration 
mechanism must treat both particle and wave dynamics on an equal footing. In fact, the situation 
is even more difficult, since the acceleration process turns out to be so efficient that the pressure 
of accelerated particles markedly modifies the structure of the shock by both the overall shock 
compression and the flow profile. Historically, these three aspects of the mechanism have been 
considered in isolation. First, a test particle theory was formulated, in which wave generation 
was only tacitly implied in the prescribed particle diffusivity, and the backreaction of accelerated 
particles on the shock structure was neglected. The latter was first included in the framework of the 
so called two-fluid model, where the accelerated particles contribute to the energy and momentum 
fluxes across the shock, but were assumed to be massless (Axford, Leer, & Skadron 1977; Drury 
& Voelk 1981). Later, kinetic models were developed, both numerically and analytically. Early 
theories that include wave generation selfconsistently with the acceleration process, in turn neglect 
particle back-reaction on the shock (Bell 1978). The importance of a self-consistent treatment of 
the nonlinear modification of the shock structure and wave propagation has been demonstrated by 
Malkov et al. (2002). That this three-way coupling (between particles, waves and flow) is indeed 
strong can be understood by considering compression of particle-generated Alfven turbulence in 
a nonlinearly modified (converging) plasma flow ahead of the shock. Since the wave number of 
the Alfven waves (which are almost frozen into the flow, Va <C Vshock) is increasing because of 
compression, particles with highest energies can no longer interact resonantly with the waves and 
so simply leave the system along the field lines. This limits the acceleration rather naturally as soon 
as particle pressure reaches a nonlinear level (comparable to the flow ram pressure), sufficient to 
modify the flow and cause a significant wave compression. Generation of longer wavelength, larger 
scale waves, as described above to initiate this process is clearly desireable. We began to study 
this process in (Diamond & Malkov 2004) and we continue this study in this paper. In the spirit of 
the discussion of this paragraph, we treat both wave and particle kinetics on an equal footing in the 
analysis presented here. 

The remainder of this paper is organized as follows. Section 2 presents the overall structure 
of the theory in the context of a discussion of accelerated particles, plasma flow and waves near 
the shock front. Section 3 discusses the dynamics of wave interactions in the shock precursor. It is 
divided into three subsections, dealing with Alfven wave turbulence, acoustic wave turbulence, and 
the dynamics of Alfvenic-acoustic coupling. Section 4 presents the theory of induced diffusion of 
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Alfven wave quanta. Section 5 deals with mechanisms of energy transfer to larger scales. Section 
6 presents conclusions and a discussion. 



2. Accelerated particles, plasma flow and waves near the shock front 

The transport and acceleration of high energy particles (CRs) near a CR modified shock is 
usually described by a diffusion-convection equation. It is convenient to use a distribution function 
f(p) normalized to p 2 dp. 

d ± + u d ±_^± = X J-IL p d ± (5) 
dt dx dx dx 3 dx dp 

Here x is directed along the shock normal which, for simplicity, is assumed to be the direction of 
the ambient magnetic field. The two quantities that control the acceleration process are the flow 
profile U (x) and the particle diffusivity K(x,p). The first one is coupled to the particle distribution 
/ through the equations of mass and momentum conservation 



Wt P + Tx PU=0 (6) 
l- tf>U + l. {pU 2 + P c+ P s)=0 (7) 



where 



P c (x) = ^mc 2 r P*P=f(p, x ) (8) 

J Pinj VP +1 

is the pressure of the CR gas and P g is the thermal gas pressure. The lower boundary pi n j in 
momentum space, which separates CRs from the thermal plasma, enters the equations through 
the magnitude of / at p = pi n j. This value specifies the injection rate of thermal plasma into the 
acceleration process. The particle momentum p is normalized to mc. Note that the two-fluid model 
can be derived from the system of eqs. (5-8) by taking the energy moment of equation (5). The 
spatial diffusivity K, induced by pitch angle scattering, prevents particle streaming away from the 
shock, and thus facilitates acceleration by ensuring the particle completes several shock crossings. 

Both two fluid and kinetic treatment of the system (5-8) indicate a marked departure from 
test particle theory. Perhaps the most striking result of the nonlinear treatment is the bifurcation 
of shock structure (in particular the shock compression ratio) in the parameter space formed by 
the injection rate, shock Mach number, and maximum particle momentum (Malkov, Diamond, & 
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Volk 2000). In particular, for sufficiently strong shocks and high particle energies the transition 
from the test particle (unmodified) shock solution to the strongly modified, efficiently accelerating 
shock solution is not gradual. It has been hypothesized that in the critical range of parameters other 
physical processes must play a crucial role. This include plasma heating and a modified particle 
confinement regime, both of which are intimately related to the wave and turbulence dynamics. 
We consider this in the next section. 



3. Dynamics of Wave Interactions in the CR Shock Precursor 

The mechanism of transfer or scattering of Alfven wave energy to larger scales is rather differ- 
ent from that associated with the conventional picture of a turbulent MHD dynamo. Most notably, 
it is a redistribution of wave energy from kr g ~ 1 to kr g < 1, where it can consequently con- 
fine higher energy particles, since kp/m = Q. at resonance. This process is not one of mean field 
generation or magnetic flux amplification, though the transfer to the large scale can generate an 
apparent B rms on those scales. Here, magnetic fluctuations are produced via the familiar process of 
cyclotron resonance of cosmic rays. The energy transfer mechanism under study is simply a decay 
or modulational instability, which is a nonlocal transfer of Alfven wave energy to larger scales, 
via scattering off density perturbations in the shock precursor. In contrast to an inverse cascade, 
this transfer is nonlocal in scale. Also, we note that the shock precursor is itself linearly unstable 
to acoustic perturbations. This mode is called the Drury instability, and assures a well populated 
scattering field of density perturbations off which Alfven waves scatter. Given the plausible as- 
sumption that the large scale scatterer field is stochastically distributed {i.e., consists of randomly 
phased acoustic waves), the effect of the decay process on the Alfven wave spectrum is to produce 
a random walk in k, via random refraction. 

The spatial structure of an efficiently accelerating shock, i.e., a shock that transforms a signif- 
icant part of its energy into accelerated particles, is very different from that of an ordinary shock. 
The most extended part of the shock structure consists of a precursor formed by the cloud of ac- 
celerated CRs diffusing ahead of the shock. If the CR diffusivity K.(p) depends linearly on particle 
momentum p (as in the Bohm diffusion case), then, at least well inside the precursor, the velocity 
profile U (x) is approximately linear in x (Malkov 1997) where x points downstream antiparallel 
to the shock normal, Fig. 2. Ahead of the shock precursor, the flow velocity tends to its upstream 
value, U\, while on the downstream side it undergoes a conventional plasma shock transition to 
its downstream value Ui (all velocities are taken in the shock frame). This extended CR precursor 
(of the size Lqr ~ ^(Pmax)/Ui) is the place where we expect turbulence to be generated by the 
CR streaming instability and where it couples to longer wavelengths. In particular, the density 
fluctuations in the precursor are what refracts the Alfven waves (produced by the same energetic 
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particles which form the precursor) and scatter them to larger scales. 



3.1. Alfven Wave Turbulence 



The growth rate of the ion-cyclotron instability is positive for the Alfven waves traveling in 
the CR streaming direction i.e., upstream, and it is negative for oppositely propagating waves. The 
wave kinetic equation for both types of waves can be written in the form 



Here (k,x,t) denote the population of quanta propagating in the upstream and downstream 
directions, respectively. Also, co^ are their Alfven wave frequencies, a> ± = kU ±Wa, where Va is 
the Alfven velocity. The linear growth (y + ) and damping (y - ) rates are nonzero only in the resonant 
part of the spectrum, for which kr g (p max ) > 1, i.e., y ± = y ± (k). In the most general case, the last 
term on the r.h.s. of equation (9) represents nonlinear interaction of different types of quanta N + 
and N~, and also compressible MHD self-interaction of each population via steepening. As seen 
from these equations, the coefficients in the wave transport part of this equation (i.e., l.h.s.) depend 
on the parameters of the medium through U and Va, which in turn, may be perturbed by slow, 
large scale fluctuations. This usually results in parametric or modulational phenomena (Sagdeev 
& Galeev 1969). We will concentrate on the acoustic type perturbations (which may be induced 
by Drury instability), so that we can write for the density p and velocity U 



p = p + p; U = U + U 
The variation of the Alfven velocity Va = Va — Vao is then 

V A ^-\v A ^-. 
2 po 

For simplicity, we assume that the plasma |3 < 1, which is valid upstream of the subshock but 
not downstream. Since we are primarily interested in the upstream, shock precursor turbulence, 
this assumption is at least reasonable, though one should address the (3 ~ 1 case as well. We will 
consider this in a future publication. In the downstream medium of a highly superalfvenic shock, 
clearly (3 3> 1. We neglect the variation of U compared to that of Va in equation (9). Note that 
fluctuations in U merely diffuse the location of the Alfven wave population in the precursor flow 
field. The perturbations of Va in turn induce perturbations of Nf, so we can write 



JV± = (N ± )+N ± 
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where (A^ (k,x,t)) is the quanitty of interest, namely the mean wave population. This is obtained 
via quasi-linear theory, applied to to the wave kinetic equation in the same way it is usually applied 
to the particle kinetic equation. Given that our goal is to obtain an evolution equation for the 
average number of Alfven quanta (A^), averaging equation (10) then yields 



l(N ± ) + (U±V i )l(N±)- <^>4 (^|^) = f + (C (10) 

Here the subscript x stands for the jc-derivatives. In order to calculate the correlator (j^N^ in 

equation (10) via quasi-linear closure, one must first determine the coherent response AP 1 (k,x,t) 
to the density perturbation field p/po- This requires solution of equation (9), i.e., the calculation 
of modulation. Since the precursor density fluctuations modulate a state of interacting, finite am- 
plitude Alfven waves, this constitutes the unperturbed (unmodulated) state. Thus, the solution of 
equation (9) is determined by a procedure similar to the Chapman-Enskog expansion. To lowest 
order, then 



YfA^+C^jA^Ar} ~0 (11) 

This relation implies that to lowest order, the linear growth rate y + is in balance with the local 
nonlinear term C + , and the linear damping of counter streaming waves y~ is in balance with their 
generation by the scattering and conversion of forward propagation waves. It is useful to note here 
that C is in general a 2 x 2 matrix operator, each component of which is nonlinear in N. Formally, 
the solution of equation (11) defines the unmodulated wave streams driven by cyclotron resonance 
of cosmic rays and inter-stream interaction, and damped by energy coupling to small scales. To 
first order, then 



^ +{u+Va) ^ _ kUx ^ _ + «^) = _f i! <w±> (12) 

gives the equation for the response of N to p x / p. The problematical element of this equation is the 
last term of the l.h.s. (y ± N ± + [8C/8N] A^). We write a 'Krook approximation' to this term as 

Y^N* + s-iV* w -AcotA^ (13) 
oA' 

where AcOk represents a nonlinear decay or damping rate. We argue for the validity of this approx- 
imation by noting that 
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a. ) since y + > 0, then C < is required for unperturbed stationarity (i.e., equation [11]) 

b. ) since C represents nonlinear interaction, C must be nonlinear in N, so N8C/8N > C. For 
example, for C = -<x/V 2 , yN - <x/V 2 = gives the unperturbed value Nq = y/cc and NdC/dN = 
-2aN N = -yN 

c. ) so, y^N^ 1 + (SC/SA^Af* < 0, and hence corresponds to a damping rate, which may be 
approximated ala' Krook as Aco^/V*. Note that this result is also consistent with the requirement 
of causality, for A(D k — > O (e). Note here that the approximate form of the matrix C is now diago- 
nalized. 

It is tempting at this stage to take Aaof ~ jy*], where the absolute value is required by con- 
sistency with causality. However, this direct balance can be established only for resonant waves, 
whereas we are primarily interested in the larger wavelength range kr g (p max ) < 1, where y^ 0. 
Thus, AcOfc is really due to nonlinear processes in that part of the spectrum. However, to the extent 
that longer waves are generated by the decay of shorter waves which are generated by resonance, 
some link between A&k and yt persists. A derivation of Aro^ which treats the effect of strong 
nonlinear refraction which occurs on kr g < 1 scales is presented later, in Sec.4. 

To calculate the refraction term we can now write equation (10), linearized with respect to 
iV^, as: 



where 

is the linear propagator with an eddy damping rate Acfl^ . Solving equation (14) for N ± , we thus 
obtain the following mean field equation for (N^) 



| <^> + 4 <^> - «4 <^> - 1*4 <^> + <c (i5) 

Notice that the mean field approximation to the refraction term in the wave kinetic equation is 
a diffusion operator in k space for the Alfven wave spectrum. This diffusion represents random 
refraction by the acoustic perturbations p via the density dependence of Va so 
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D k is an example of the well-known phenomenon of induced diffusion. Induced diffusion is a 
generic type of 3-wave interaction due to resonant triads with k + k'+q = but |k| , |k'| 3> |q|. 
Such triads can be represented by thin, isosceles triangles, Fig. 3. The basic physics of induced 
diffusion is random refraction by large scale perturbations. Induced diffusion can thus be obtained 
from eikonal theory approaches. The inversibility which ultimately underpins the diffusion has its 
origins in the stochasticity of the Alfven wave rays in the field of density perturbations. Trans- 
forming to Fourier space, we first represent p as 



p = £p/^-*V 

q 

and note that due to the local Galilean invariance of L, we can calculate its Fourier representation 
in the reference frame moving with the plasma at the speed U (x) as: 



L 



£ A = ±iqV A + to£-kU x l (17) 



Then, equation (16) can be re- written as: 



Pq 

Po 



(18) 



The last (wave refraction) term on the r.h.s. of equation (17) can be estimated as Uf /K(p max ), 
which is the inverse acceleration time and can be neglected in comparison to the frequencies qV A 



and A(0. Hence, for 9tL, „ we have 

K,q 



Acq 



(19) 



Note that equation (19) states that the correlation time x c for induced diffusion is set by the Alfven 
wave damping time, (~ l/Acfl^) reduced by the effects of finite transit time for a wave to propagate 
through the density perturbation, if % c qV A > 1. Clearly, equation (19) has two limits, a 'strong 



turbulence' limit where 9tL 

,2 T/ 2„ 



±-1 
'k,q 



AG)" 



T cy fc and a 'weak turbulence limit' where 9tL 



±-i 

k,q 



l/q V^x c k. For further convenience, we write the population density of acoustic waves (phonons) 
as 



" o>' q 

where W q is the energy density of acoustic waves (with to* = qC s their frequency). 
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W -C 2 ^ 
Po 



For Dk in equation (15) we thus finally have 

fl*=|^i>X 2 jf\ ±2 % (2°) 

Note that represents the rate at which the wave vector of the Alfven wave random walks due 
to stochastic refraction. Of course, such a random walk necessarily must generate larger scales 
(smaller k), thus in turn allowing the confinement of higher energy particles to the shock. Thus, 
confinement of higher energy particles is a natural consequence of random Alfven wave refraction 
in acoustic perturbations. 



3.2. Acoustic turbulence 

In contrast to the Alfvenic turbulence that originates in the shock precursor due to cyclotron 
emission from accelerated particles, there are (at least) two separate sources of long wave acoustic 
perturbations. One is related to parametric and modulational (Sagdeev & Galeev 1969; Skilling 
1975) processes undergone by the Alfven waves. These take the usual form of decay of an Alfven 
wave into another Alfven wave and an acoustic wave. The other source is the pressure gradient of 
CRs, which directly drives linear instability. The latter leads to emission of unstable sound waves 
via the Drury instability. By analogy with equation (10), we can then write the following wave 
kinetic equation for the acoustic waves: 

3 3 3 

-N q + U Y N q - q U xTk N q =(f q + y D )N q + C{N q } 

Here Jd is the Drury instability growth rate Jo = Jo [VPcr] and 7^ is the growth rate of the decay 
instability r fq= n f q l [N^]. We first consider the decay instability of Alfven waves. Note, however, 
that the combination of Drury instability and decay instability can lead to more rapid generation 
of mesoscale fields by coupling together the linear and nonlinear processes, both of which amplify 
density perturbations. 
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3.2.1. Modulational Instability ofAlfven Wave Packets in a Density Scatterer Field 

The mechanism of this instability is the growth of the density (acoustic) perturbations due 
to the action of the ponderomotive force on the acoustic waves by the Alfven waves. This force 
can be regarded as an effective radiation pressure term which must appear in the hydrodynamic 
equation of motion for the sound waves (written below in the comoving plasma frame) 

Eliminating velocity by making use of continuity equation 

dp dV n 

we obtain 



The Alfven wave pressure can be related to their population densities via the energy density, so 

k V 7 

Using the relation (14) between the density perturbations and the Alfven waves and separating 
the populations of forward and backward propagating sound waves p*, we can obtain from equa- 
tion (21) the following dispersion relation for the nonlinear growth of the the acoustic branch via 
modulational interactions 

» 2 -? 2 c? = 9 2 I^vii^'|<^> 

upon writing G) = ±qc s + iy ± , we then find the following growth rate of acoustic perturbations 

Note that modulational instability requires an inverted population of Alfven quanta ie d (N) /dk> 
0. As Alfven waves are generated by high energy resonant particles in a limited band of k space 
at short wavelength (i.e., kr g ~ 1), such an inversion clearly can occur. Note also that shorter 
wavelength modulations appear to grow faster, but this trend saturates when q > Ago^/V^. Note 
also that the coherence time between the Alfvenic packet and the modulating density perturbation 
field is set by L k l q ~ AtOjt/ (q 2 Vg r + A(xfy which, not surprisingly, also sets the correlation time in 
D k . 
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3.2.2. Drury Instability 



There is a linear instability which also leads to efficient generation of acoustic waves and is 
driven by the pressure gradient of the CRs in the shock precursor. Such a VPc^-driven process 
is of particular interest, as it taps free energy which is generated and stored as a consequence of 
the acceleration process, itself. The growth rate has been calculated by Drury (1984) and Drury & 
Falle (1986) (see also Zank, Axford, & McKenzie, 1990; Kang, Jones, & Ryu, 1992), and can be 
written in the form: 



Here Pq and Pq x are the CR pressure and its derivative, respectively, and jc is their adiabatic 
index. Note that instability requires dhiK/dlnp > —1, so that the structure and dependencies of 
the cosmic ray diffusivity are critical to the Drury instability. The physics of this instability is 
that the density fluctuations in the decelerating precursor flow induce the CR pressure gradient 
fluctuations which (unless they are exactly proportional to the density fluctuations) make the flow 
deceleration non-uniform and so can amplify the initial density fluctuations ( i.e., Drury, 1984; 
Kang et al., 1992). For an efficiently accelerating shock, the adiabatic index jc in equation (22) is 
y c m 4/3. Note that we have omitted a term — U x which is related to simple compression of wave 
number density by the flow, and which would enter the r.h.s. of equation (21) (see Drury & Falle, 
1986). This term is smaller, by a factor of C s /U , than the second (destabilizing) term. The first 
term is damping caused by CR diffusion, calculated earlier by Ptuskin (1981). 



Compressibility of the media in a shock environment is believed to be responsible for the 
formation of coherent magnetic structures observed upstream of the Earth's bow shock and in- 
terplanetary shocks. These structures are thought to evolve from phase steepened Alfven wave 
trains and should be relevant to the SNR shocks as well. The bottom line of this highly evolved 
and mature field of research is that parallel compressibility transforms Alfven wave trains into 
steepened Alfven solitons, described by the DNLS (Derivative Nonlinear Schrodinger Equation) 
(Kennel et al. 1988a; Medvedev 1999)and its variants. The DNLS describes, within the framework 
of coherent interactions, the refraction of Alfven waves by density perturbations and its feedback 
on the Alfven wave envelope. Thus, the DNLS describes a process which is the coherent analogue 
of k- space diffusion and which does not suffer from the limitations of eikonal theory. A simpli- 
fied derivation of the DNLS is presented below. Consider the familiar Alfven wave dispersion law 
(0 = &||Va- It stems from the differential equation for the magnetic field perturbations. If we now 




(22) 



3.3. Nonlinear Wave Trains and Alfven-Acoustic Coupling 
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entertain the possibility of slowly varying density perturbations, we have 

dt 2dz A po ^ ^ 

where B is the magnetic field perturbation envelope function in the Alfven wave, and p is the 
density perturbation induced by parallel compression. From the continuity and parallel momentum 
equations we have 



¥ = - p °!F <24) 

*n£ = -4» + 8i+V) <25) 

with pV 2 /2 ~ B 2 /8n (i.e., both V and B enter ala' Bernoulli) and V ~ B since the basic waves 
are Alfven waves. Assuming that the perturbations propagate at approximately the Alfven veloc- 
ity, from equations (24-25) we obtain the relation between wave envelope perturbation levels and 
density fluctuation levels, i.e., 

Po 2 V B 2 

Note that this relation gives a clear and simple relation between density and magnetic fluctuation 
energy levels. Substituting this in equation (23) we obtain 



9B 13/ V A 
+ 



B 



B = 



dt ' Adz \ \-C 2 /Vl) B 2 

Dispersive corrections to the Alfven mode can be easily added into this equation as well. These 
ultimately limit steepening, producing collisionless shocks. For the singular case |3 1, a modi- 
fied kinetic equation (KNLS) can be derived (Medvedev & Diamond 1996). The DNLS equation 
possesses soliton solutions that may self-organize in quasi periodic structures (wavetrains). Upon 
addition of dissipative term and a driver to this equation (e.g., growth due to instability), quasi 
periodic shock train solutions may also be obtained(Kennel et al. 1988b; Hada et al. 1990). Such 
DNLS and KNLS shocks can trap and mirror energetic particles, thus enhancing their confine- 
ment. Consideration of DNLS/KNLS structures and their dynamics gives a clear and compelling 
physical picture of the Alfvenic mirroring process as well as constituting a coherent analogue of 
the stochastic modulation discussed above. In particular, the coherent modulation of instability 
amplifies the magnetic energy of the envelope, and thus couples energy to large scales. 



- 19- 



4. Wave Quanta Evolution in Strongly Refracting Turbulence 

Noting that Alfven waves are generated at kr g ~ 1 , while the modulational interaction which 
scatters the spectral population to large scales occurs at kr g < 1 . In the likely case that this inter- 
action is strong (i.e., p/po ~ 1) and stochastic, it can directly produce wave packet decorrelation, 
which also contributes to the interaction damping rate A(0 fc . Here, we calculate this decorrelation. 
Since in this case, stochastic refraction is dominant, the details of the l.h.s. of the wave kinetic 
equation are not important, so we write, in the precursor flow frame 

dN± , . . s dN± d (u) dN± 3 3 p x ~ + 

IT + ««> ±v ") "aT ~3x ~aF + a;"" + lx kv *t N = (26) 

The mean field equation is just 



Px 



(27) 



Note that the problem of determining the evolution of the mean wave populations reduces to cal- 
culating the correlations (uN^ and ^(p^/po) N^, which constitute the spatial and wavenumber 

fluxes, respectively. In the spirit of quasilinear theory, we calculate the value of /V ± to close the 
correlators by determining the response of A^ ± to the acoustic wave perturbations u and p x /po- 
Writing the linear propagator 



(I \±\T\ 3 f 3 (") 3 T 



(28) 



the response is given by 



+ L±N±+C 



± 



= — u- 



a(/v ± ) kv APx d 



dx 



2 po dk 



Note that the Hamiltonian structure of eikonal theory means that the flow in (x, k) space is incom- 
pressible, thus allowing re- arrangement of the order of the derivatives. Here C 'V ± is that portion 
of the nonlinear terms which is phase coherent with the population fluctuation. Taking modulations 
to have the form 



( u \ 
P/P± 



( 



U q \ 



P?/po 



i(qx—Q.t) 
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then the response equation becomes 



-i(n- qLq )N^+c 



N 



± 



~ kV A . p q 3 , 

= —u a — —iq——(N) 

q,Q q dx 2 p dk x 1 



Here we have re- written L as iqL q for the corresponding Fourier mode. C \N ± is given ex- 

L J q.Q, 

plicitely by (we introduce a "two-vector" q = (q, CI) for short) 



l(uN+) + l- k kV A 



= E 



-iqu^m.- + ——Uf—±.tf 



q ' dk 2 ^ po 



q'+q 



The expression for C 



± 



may, in turn, be obtained as in quasilinear theory by approximating 



Ng +d in terms of the beats at q + q' determined by the wave kinetic equation, i.e., 



-i [Cl + Cl'- {qL q + q'L q ,)] /V± + Ag)^/V± = iqu^Nf + ^'y^f 
A formal solution then yields 



N^f = R q+q -> [iqujN? + -^iq'^^N, 



± , kV A .jp^d_~ ± 
Podk * 



where 



R-lj = -i [fl + Cl'- (qL q + q'L q ,)] + Aw q+q , 
This in turn gives the result 



n - 



Aa q Nf 

£ (-iqU-f - ^Y^V-q* J R q+q< 



Ignoring cross terms the expression for A&qN^ may be simplified to 



d dN„ 
± _ at n 1 



where 



D q = H\ U q'\ ^ R q'+q 

q' 




+q 



correspond to wave packet diffusion in position and k, respectively. Note also that 



[n + n'-(qL q + q'L q ,)Y+A(d 2 - +ir 



so that A& q is defined recursively, as usual, and that both types of scattering contribute to the total 
decorrelation. Finally, again ignoring cross terms, we can now obtain the mean field equation, 
which is 



a( M )3(^v ± ) 



dt 



dx 



-k 



d_ 

dx 



d_ 

dx 



dx 



dk 







Here D and again correspond to spatial and wavenumber diffusion of an Alfven wave 
packet, and are given by lim^o A^Afc,^ respectively. The output of this relatively straightfor- 
ward calculation is the set of renormalized quasilinear diffusion coefficients D, D^. In each of 
these the decorrelation rate is set by decorrelation rate A(i) q , since the dissipation in the acoustic 
wave spectrum is weak. Since Aco is defined recursively, explicit relations between Ag)^ and the 
refracting field spectrum are best discussed in various simple limits. In the case that decorrelation 
is dominated by spacial scattering, 
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Here Q. absorbs the propagator Doppler shift. In this case then, Acq ~ (q 2 (S 2 ) - Q 2 ) 1/2 . Note 
that a critical level of (u 2 ) is required for irreversibility, i.e., to make a stochastic Doppler shift 
sufficient to overcome the frequency Cl. In the limit where wave-number scattering dominates but 
the waves are weakly dispersive (i.e., dV g /dk — > 0) 



k 2 V 2 
K V A -2 



(AkY 



Po 



Aco„ 



A 2 + Aco 2 



so 



Aco 



k 2 V 2 



(Vp)' 



PS 1 



a 2 



1/2 



Again, a critical rms density gradient fluctuation level is required. In the limit of stronger dispersion 
and strong turbulence, the decorrelation rate is 



Aco ? ~ 



q% 




2„,ik 2 V 2 



1/4 



Note that in this case Aco ~ (Vp/po)^. It is important to stress that here, the physical processes 
represented by A& q are random advection and refraction of wave packets by the acoustic wave fluc- 
tuations on scales q, where qr g <C 1 , and not the interaction of Alfven waves with kr g ~ 1 . Since, 
of course, both processes occur, and since the modulational interaction of Alfven wave packets 
with density perturbations generates larger scale waves, the physics of wave packet decorrelation 
is, not surprisingly, strongly dependent on scale. While this may result in some technical difficul- 
ties in qualitative calculations, it does ensure that irreversible 'modulational turbulence' dynamics 
persists over a broad range of scales, and is not sharply localized at kr g ~ 1 . 



5. Mechanisms of Magnetic Energy Transfer to Larger Scales 

As should be clear from the considerations discussed above, there are a variety of nonlinear 
processes that can lead to the transfer of magnetic energy (generated by accelerated particles in 
form of the resonant Alfven waves) to longer scales. First, as it can be seen from equation (15) 
(i.e., last term on the l.h.s.), scattering of the Alfven waves in k space due to acoustic perturbations 
transfers magnetic fluctuation energy away from the resonant excitation region to smaller (and also 
to larger) k, and also amplifies the long wavelength acoustic scattering field. Second, the nonlinear 
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interaction of Alfven waves and magnetosonic waves represented by the wave collision term on 
the r.h.s. can drive such a process. Next, as we discussed in the last subsection, solutions of 
DNLS-KNLS equations can be well represented by quasiperiodic wave structures. The interaction 
between such structures, such as, for example, shock waves in the Burgers model, leads to their 
coalescence, which in turn, means an efficient transfer of excitation to larger scales. Finally, we 
remind the reader that even within the frame work of weak turbulence theory, induced scattering 
(i. e., nonlinear Landau damping) of Alfven waves on thermal protons leads to a systematic decrease 
in the energy of quanta which, given the dispersion law, again means energy scattering to longer 
wavelength. 

Continuing our main discussion of the wave refraction by acoustic perturbations generated by 
the Drury and modulational instabilities, it is important to emphasize the following point. As seen 
from the instability growth rate, equation (22), the Drury instability growth rate is proportional 
to the gradient of P c . While one would naively expect Drury instability to relax the gradient that 
drives it, we note that the outcome of the dynamic feedback we discuss here suggests that some 
of the relaxation would be off-set enhanced confinement and acceleration. Indeed, this might 
ultimately reinforce the instability, possibly triggering bursts or cyclic growth. This might help 
realize mechanisms of regulation of P c suggested earlier in (Malkov et al. 2000). 

We proceed, however, with a simpler approach by not treating the acoustic perturbations fully 
self-consistently i.e., we do not consider the connection between acoustics and particle acceleration 
and shock modification. Instead, at this level we treat them as a developed to some quasi- stationary 
state as a result of e.g., Drury instability and we consider how the Alfvenic turbulence then evolves 
in their field. Thus, we address only the simplest and most fundamental problem, that of the 
evolution of the Alfven wave population in a field of ambient density perturbations. 

5.1. Estimates of Alfven Wave Diffusion in k space 

The spectral evolution of the Alfven waves generated by accelerated particles is crucial for 
confinement and further acceleration of these particles. It is described by equation (15) and in- 
volves several processes: 

1 . wave generation at the rate y 

2. convection of the waves to the shock from the upstream side at the flow speed U 

3. blue-shift of the waves to short scales by the flow compression (term U x , see Malkov et al., 
2002) 
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4. nonlinear interaction (term C (N)) with each other and particles (induced scattering) 

5. random scattering in wave number on acoustic perturbations (diffusion in k) 



Clearly, an exact solution of the problem is very difficult, particularly because all these phenomena 
are coupled and related and, ideally, should be treated self-consistently. For example, the wave 
turbulence level N(k) directly enters K(p) of the convection-diffusion equation (5), the solution 
of which determines both the flow profile U (x) (through P c and eqs.[6-7]), and the growth rate y 
(through P c ) and also the Alfven quanta scattering rate D k (through the Drury or parametric insta- 
bilities). Each of these appears explicitly in the wave kinetic equation, thus closing the feedback 
loop. A fully self-consistent level of description would jump several steps ahead of the the cur- 
rent DSA models and is clearly beyond the scope of this paper. Therefore, to achieve any simple 
understanding, several simplifications are necessary. In particular, we wish to focus here on the 
spreading of the wave population in k space, assuming that the spectrum is initially generated by 
accelerated particles upstream of the shock in the resonant k domain. Therefore, we ignore nonlin- 
ear interaction between waves, as well as wave refraction described by the flow compression U x . 
This effect, and (partially) also the wave self-nonlinearity have been considered by Malkov et al. 
(2002). 

First, we rewrite equation (15) in the following simplified form 



^ + U — — —D — — — — (29) 
dt dx dk k dk Wa dx 

where we have introduced the dimensionless wave intensity 4 normalized to the background mag- 
netic energy 



\B k \ 2 



along with the partial particle pressure P(x,p) in the p 3> 1 range (see equation [8] ) normalized 
to the upstream ram pressure piWj 



p tomcat (3Q) 

We do not distinguish here between the forward and backward backward propagating waves or 
Alfven and magnetosonic waves. Rather, we simply incorporate all types of magnetic fluctuations 
in the spectral density 4- The wave induced diffusivity from equation (20) can be represented 
as 
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^ = 7^aT J 2 ^ a / 4 (3D 

To further simplify this expression we note that there are several reasons for the resonance broad- 
ening A(0 in the last expression, as discussed earlier. One is the nonstationarity of the Alfven waves 
being scattered by the acoustic fluctuations, which yieldsAd) ~ y. A different reason is the refrac- 
tive scattering itself. It comes in two flavors, as discussed earlier, one of which is related to the 
dispersion of the group velocity, V g = dV g /dk, acting in concert with the scattering in k, namely 

A(0 ~ (qV'g) 2 ^ ' D 1 ^ 3 ', and is analogous to the resonance broadening familiar from wave-particle in- 
teractions in plasmas (Dupree 1972). This effect is small for nearly parallel propagation of Alfven 
waves, since V g const . However, for a more realistic spectrum composed of a mixture of Alfven 
and magnetosonic waves with a sizable off-axis component, the effect is significant. The second 
flavor is simply direct scattering of the Alfven quanta in k, regardless of V g variations. This con- 
tributes to the resonance broadening as A(0 ~ D^/Ak 2 , where Ak is the width of the wave packet. 
The first effect is usually larger in most situations when the field of scatterers is weak (p q <C po) 
but in our case, rather the opposite is true (Kang et al. 1992), so that the second effect turns out to 
be stronger. This is ultimately due to the weak dispersion of quasiparallel Alfven waves. 

For steepened acoustic fluctuations, we can assume the scatterer field is simply an ensemble 
of shock-like discontinuities, so 




This form is appropriate for a shock ensemble with a characteristic amplitude ~ y/~A ~ Ap/po 
and a separation ~ \/qo- This kind of density perturbation structure is to be expected from the 
nonlinearly developed Drury instability (Kang et al. 1992) since that would produce a pattern of 
discontinuities. From equation (31) we then obtain 



where y 3 = DkV' 2 / qoV g . Further progress is possible only for the limiting cases of weak and strong 
scattering. Therefore, in the case of strong scattering (y> 1) we have 



371 V} n o 

D k ^—Aq ^-k 2 = Vk 2 (34) 
4 V g 

where V' cancelled out. In the opposite case of weak scattering (y <C 1) we obtain 
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Restricting consideration to the case of strong scattering with as given by equation (34), the 
steady state limit of equation (29) becomes 



—k 2 — - — — - - 2M A — — 
dk dk V dx Vk dx 



(36) 

pk—m(S) c 

In the case of weak diffusion, with no spectral spreading (i.e., v — > 0), the turbulence level is simply 
directly proportional to the particle partial pressure 



/ ~ 2M A k~ l P (p = & c /kc) (37) 
where p is normalized to mc. It is convenient to introduce the following dimensionless variables 

x=— , Q = 2M A P/k, v = \n(kc/(D c ) *¥ = Ie% + * 
so that the equation for *P becomes 

d*¥ d 2x ¥ v,x3(2 

= e 2 + 4 — 

dx dv 2 dx 
This can be solved to obtain the wave spectral density / 

J-°°y/4%(%-T>) i-oo 3X /V ^ 

By introducing the dimensionless coordinate £ = x/L, it easy to see that the main parameter that 
determines the behaviour of the solution for I is the stregth parameter S 

g vL t com 



s 



where t conv = L/u\ is the precursor (of length L) crossing time and t s = 1/v is the refractive scat- 
tering time for the wave. In the limit 5^0 one simply obtains I = Q, equation (38), which 
is equivalent to equation (37) above. In the more interesting case S 3> 1, by using the steepest 
descent method we obtain 

'^L^-^^left'.v')^.,,,.,, (39) 
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where H is the Heavyside function. One sees that there is a S 3> 1 reduction in the wave spectral 
density as compared to the scattering free case given by equation (37). To clarify this last result, let 
us specify the particle distribution function at some dimensionless distance £, ahead of the shock. 
We assume that the upstream medium is at x < and since we neglected the flow compression in 
equation (36) we assume the stationary solution of the diffusion-convection equation, equation (5) 
has the following form 



/(p,x)=/ (p)exp(^) (40) 

Note that in the case of a modified shock precursor, the solution has a similar form, except for the 
modified term in the exponent. That is replaced by an integral over x and a multiplicative factor of 
the order of unity, depending slowly on p (Malkov 1997). For this simple estimate we can assume 
that K is a linear function of p, i.e., 



K~K fc — (41) 
Pb 

where we have introduced the lower cut-off momentum pb from equation (40) in an obvious man- 
ner, K/, (x) = K[pt, (x)] = \x\ u\. We also assume that the spectrum has an upper cut-off at p max and 
K scales linearly with momentum up to p max , so that K max ~ Kb (p max / pb). We account for the fact 

thatL = Kmax/Ui. 

From the point of view of confinement and thus acceleration improvement due to spectral 
transfer it is important to understand whether the spectrum that is resonantly generated by a group 
of already accelerated particles near p max propagates to lower k, so as to facilitate confinement and 
acceleration of particles with p > p max - Note that the wave compression due to the nonlinear shock 
modification in the CR precursor leads to an opposite evolution in k. 

Since the waves are originally driven by accelerated particles, it is instructive to express the 
wave energy density through the particle pressure. Note that there is a coordinate dependent low 
energy cut-off pb (x) on the particle spectrum which physically means that low-energy particles 
cannot diffuse far ahead of the subshock (see equation [40]). Normalizing the particle pressure to 
the shock ram pressure it is convenient to express the former one as an acceleration efficiency 



e (x) = = I . . dppifo ( P ) (42) 



P c (x) 4k mc z fP max A 3 
piwf 3 piu\Jp b ( x ) 

where we have ignored the exponential factor in the distribution function, equation (40) since the 
integral in p effectively runs above pb{x). Assuming then that pb is reasonably close to Pmwc, 
Pb ^ Pma X , which means that x is taken sufficiently far from the sub-shock, but not so far so there 
are not enough particles or waves. This case is therefore interesting, in that there is still sufficient 
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space for spectrum spreading before waves are convected into the subshock. In this limit, we can 
estimate the spectrum from equation (39), expressing it for convenience in terms of the resonant 
p, rather than in terms of k, as: 



, W = ^e W ( /' (43) 

P< Pmax ~ Pb 

The connection of this result to scattering of the Alfven quanta in k can be seen directly from 
equation (36). It corresponds to the solution of equation (36) in two regions of wave number 
k, where the source is absent, P ~ 0. The low-fc asymptotics (high p) is simply a solution with 
constant / (k) = const , whereas the high-A: part of the solution corresponds to the constant flux 
of the wave density in k space, k 2 dl/dk = const . These limiting forms obviously match in the 
region where waves are generated, i.e., where p ~ p b ~ p max - Obviously, equation (43) indicates 
that diffusive refraction in k space leads to the confinement and acceleration of higher energy 
particles, with p > p ma x- 

At this point, it is appropriate to discuss, more generally, the flow of wave population density 
N and wave energy density £ in k space which results from the stochastic scattering of Alfven 
waves off precursor density fluctuations which we described above. Since this interaction is mod- 
eled (at the quasilinear level) by diffusion in k space, we can write the relevant part of the wave 
population kinetic equation as 

d(N)_ d d(N) 
dt dk k dk 

so a gradient in the population (N) produces a flux in k 

a <A0 



n = -D k - 



dk 



Thus, since waves are excited at kr g ~ 1, d (N) /dk > for kr g < 1 and d (N) /dk < for kr g ^> 1 . 
Thus, < for kr g < 1, suggestive of a flux of the Alfven wave population density to large scales. 
The total wave population 



niax 



91 = I dk(N) 

evolves according to 



k>nax 
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and thus is determined by the slope of the spectrum N in the high and low k limits since the wave 
energy density •£ = CO^/V, the evolution of the total energy 



A detailed spectral evolution calculation, beyond the scope of this paper, is required to precisely 
determine the sign of dE / dt. That said, it seems very likely indeed that the wave number integrated 
contribution (second term on r.h.s.) is both dominant, and is negative. This follows from > 0, 
V g = doo/dk > (i.e., the counter-propagating streams do not precisely balance) and d (N) /dk > 
0, since there is a greater population density at small scales. Thus dE/dt < 0, which is also 
consistent with the outcome of the modulational instability calculation. Since waves are scattered 
by refraction by precursor density fluctuation, we know that the nonlinear interaction leaves 



precursor density perturbation energy and depletes Alfven wave energy while it scatters Alfven 
wave energy towards larger scales. 

The upshot of this decay/modulational interaction process can thus be viewed as an energy 
transfer process in k space which increases the population of high energy cosmic rays which are 
confined to the shock and precursor. This may be seen by considering the sequence below, i.e., 

i. ) as usual, energetic CRs generate Alfven waves with kr g ~ 1 

ii. ) these waves scatter off of ambient density perturbations in the precursor via decay insta- 
bility, thus: 

a.) producing larger wavelength Alfven waves 




is given by 




4(£w+£ P )=0 



so dEp/dt > 0, i.e., the decay instability of short wavelength (kr g ~ 1) Alfven waves amplifies the 



b.) amplifying the density perturbation field and so producing a flow of fluctuation energy 
(including magnetic energy) toward longer scales. At the same time: 
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a. ) the longer wavelength Alfven waves will confine higher energy particles 

b. ) the precursor density perturbations will produce high efficiency acceleration (i.e., beyond 
the Bohm limit) of particles with p > p% (where p.* corresponds to the "knee") by scattering these 
particles in a converging flow. This process does not require particles to cross and re-cross the 
shock, itself. 

In this way, we see a direct connection between fluctuation energy flow to larger scales and 
achieving increased particles energy during the acceleration process. Note also that two accel- 
eration processes {i.e., the traditional one involving multiple shock crossing via confinement by 
turbulent pitch angle scattering and the recently proposed one involving scattering off inhomo- 
geneities in the converging percursor flow) can work together by exploiting the wave scattering 
mechanism discussed here. 



6. Conclusions and Discussion 

In this paper, we have examined the dynamics and generation of mesoscale magnetic field in 
diffusive shock acceleration. The principal results of this paper are: 

1. the identification of a new (in the context of DSA theory) nonlinear, multiscale interaction 
mechanism involving precursor density fluctuations and resonantly generated Alfven waves 
was identified. This mechanism is one of modulational or decay instability of Alfven waves 
in a scattering field of density perturbations. 

2. the demonstration that such interaction: 

(a) generates larger scale waves, which in turn confine higher energy particles to the shock, 
thus allowing their acceleration 

(b) can also enhance the level of density fluctuations in the precursor, and so assist in the 
(recently proposed) high efficiency acceleration by scattering cosmic ray particles off 
inhomogeneities in the converging flow. 

3. the explicit construction of a strong turbulence theory, based on modulational interaction, 
for calculating the spreading of Alfven wave energy. The principal mechanism by which 
Alfven wave energy is scattered to larger scales is shown to be random refraction by the 

spectrum (Vp/p) - of density gradient fluctuations in the precursor flow. These density 
fluctuations can be produced by the modulational process itself, or develop from precursor 
instabilities, such as the Drury instability. 
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4. the identification of the critical parameter which governs this k space diffusion process, 
namely S = L/u\X s , where \/% s ~ D^/Ak 2 is the refractive scattering time. 

5. the demonstration that for S > 1, diffusion in k will generally broaden and suppress narrow 
band spectra at kr g ~ 1, while scattering energy to longer (and shorter) wavelengths. 

The relation of this type of mechanism to existing concepts in DSA theory was discussed thor- 
oughly. 

The results of this paper have several interesting implications for high energy cosmic ray ac- 
celeration. First, given that density fluctuations are always present via the Drury instability, waves 
at kr g ~ 1 will not grow to large amplitude (85 > Bq), but rather will have their energy diffused in k 
to a broad band of larger and smaller scales. Note that this new nonlinear mechanism of saturation 
of the resonant RMS magnetic energy is independent of the quasilinear-type instability saturation 
mechanisms discussed in the Introduction, which were considered earlier by e.g., MacKenzie & 
Voelk (1982) and Achterberg & Blandford (1986). Thus, theories which predict the generation of 
strong, small scale fields without considering scattering by precursor fluctuations probably have 
significantly overestimated the strength of the field at kr g ~ 1 . Second, the modulational mech- 
anism presented here, while not a "dynamo", in a strict sense, is a robust and universal means 
to scatter wave energy to larger scales and so confine higher energy particles {i.e., p > p m ax)- 
Hence, it constitutes a novel means for enhanced acceleration. A quantitative calculation of the 
resulting energy spectrum clearly requires solution of the coupled kinetic equations for Alfven 
waves, acoustic perturbations, along with the evolution equation for the energetic particle distri- 
bution function /, and the nonlinear shock conditions. This is obviously a formidable task, and 
one which will require significant effort in the future. However, a more tractable approach i.e., a 
first step is to solve only the coupled equations for the resonantly generated Alfven wave popu- 
lation and the energetic particle distribution function, incorporating diffusion in wavenumber by 
a prescribed spectrum of ambient density fluctuations associated with Drury instability. This cal- 
culation, which is a fairly modest extension of the analysis in Sec. 5.1 of this paper, should reveal 
the quantitative relations between the precursor scatterer field spectrum and the energy spectrum 
of accelerated particles confined to the shock by turbulent pitch angle scattering. The results of 
this calculation will be presented in a future publication. However, the reader should also keep in 
mind that the modulational instability will amplify any perturbations in the precursor, thereby also 
enhancing acceleration by scattering of CRs by inhomogeneities in its converging flow. Thus, it 
increasingly seems that the most energetic particles result as much from precursor dynamics, as 
from the traditionally invoked crossing of the subshock discontinuity, per se. 

We thank L. Drury, H. Li and M. Medvedev for stimulating discussions. This work was 
supported by NASA under grant ATP03-0059-0034 and by the U.S. DOE under Grant No. FG03- 
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Fig. 2. — The structure of a CR modified shock, with the flow profile U (. 
distribution P c (x). The sub-shock is at x = and the CR precursor of the 
upstream, x < 0. 



x) and CR pressure 
length L p is formed 



■q 



Fig. 3. — Three-wave interaction of two nearly opposite Alfven quanta k and k 7 with a long wave 
acoustic wave q. 
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Fig. 4. — Formation of inverted population of Alfven quanta under the condition of localized (at 
kr g ~ 1) driver. 



